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ABSTRACT 

We present the results of X-ray observations of the well-studied Te V blazar Mrk 42 1 with the Suzaku satellite 
in 2006 April 28 . During the observation, Mrk 42 1 was undergoing a large flare and the X-ray flux was variable, 
decreasing by - 50 %, from 7.8 x 10"'" to 3.7 x 10"'° erg s 'cm ^ in about 6 hours, followed by an increase 
by 35 %. Thanks to the broad bandpass coupled with high-sensitivity of Suzaku, we measured the evolution 
of the spectrum over the 0.4-60 keV band in data segments as short as ^ 1 ksec. The data show deviations 
from a simple power law model, but also a clear spectral variability. The time-resolved spectra are fitted by 
a synchrotron model, where the observed spectrum is due to a exponentially cutoff power law distribution of 
electrons radiating in uniform magnetic field; this model is preferred over a broken power law. As another 
scenario, we separate the spectrum into "steady" and "variable" components by subtracting the spectrum in the 
lowest-flux period from those of other data segments. In this context, the difference ("variable") spectra are 
all well described by a broken power law model with photon index F ^ 1.6, breaking at energy etrk — 3 ke V 
to another photon index F 2.1 above the break energy, differing from each other only by normalization, 
while the spectrum of the "steady" component is best described by the synchrotron model. We suggest the 
rapidly variable component is due to relatively localized shock (Fermi I) acceleration, while the slowly variable 
("steady") component is due to the superposition of shocks located at larger distance along the jet, or due to 
other acceleration process, such as the stochastic acceleration on magnetic turbulence (Fermi II) in the more 
extended region. 

Subject headings: acceleration of particles — BL Lacertae objects: individual (Mrk 421) — galaxies: jets — 
X-rays: galaxies 



1. INTRODUCTION 

Blazars, a sub-category of Active Galactic Nuclei (AGN), 
are characterized by broadband non-thermal emission and vi- 
olent variability on time scales from a fraction of an hour to 
a few days with strong spectral evolution. This behavior is 
best described as their emission arising from Doppler-boosted 
relativistic jets which in turn dom inates over the thermal 

] signatures seen i n other AGN (e.g., lUrrv & Padovanilll995l; 

• lUlrich et"anil997D . The broadband spectra of blazars consist 
of two peaks, one in the radio to optical-UV range (and in 
some cases, reaching to the X-ray band), and the other in the 
hard X-ray to 7-ray region. The high polarization of the ra- 
dio to optical emission suggests that the lower energy peak is 
produced via the s ynchrotron process by rela tivistic electrons 
in the jet (see, e.g.. lAngel & Stockmanlll980h . The higher en- 
ergy peak is believed to be due to Compton up-scattering of 
seed photons by the relativistic electrons. A number of blazars 
with the synchrotron peak located in the X-ray range have also 
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been detected in the TeV band, and the recent advances of TeV 
Cherenkov telescopes have revealed the existence of signifi- 
cant n umber of such TeV emitting blazars (see IWagner et alj 
I2OO8I for recent synoptic study). These objects provide the 
most direct evidence of efficient particle acceleration up to 
TeV energies, potentially providing the key information about 
the jet's composition and power, and thus the connection of it 
the central engine. 

Mrk 421 is one of the nearest (z=0.031) and brightest 
TeV 7-ray emitting blazars: it was fo und as the first ex - 
tragalactic TeV 7-ray emitting object dPunch et al.l Il992h . 
and has been repeatedly confirmed a s a TeV source by var- 
ious ground-based telescopes (e.g., lAharonian et alJ 120051 : 
Albert et al. 2007). It has also been one of the most ex- 
tensively studied blazars, and has been a target of sev- 
eral multi- wavelength campaig ns (e.g.,|T akahashi et al .ll 19961 
2000; Tost i et al.1 [1998; Rebillotetair2006; FossadeTaD 
l2008i: iLichti et al.l 1200 8). Detailed studies during large am- 
plitude flares are particularly important as they give us vital 
information allowing studies of the underlying physics of the 
source. In particular, the simultaneous observations in X-ray 
and TeV 7-ray band revealed a strong correlation between 
them and depict the picture in which the non-thermal distribu- 
tion of relativistic electrons accelerated up to TeV energies are 
responsible for the variability of both bands. Such a scenario - 
where lower energy component is synchrotron radiation from 
the relativistically accelerated electrons and the higher energy 
one is due to Compton up-scattering of the synchrotron pho- 
tons by the same electron populations themselves - is widely 
adopted as the most plausible emission mechanism for TeV 
blazars, and is cafled one-zone Synchrotron Self-Compto n 
(SSC) model (llnoue & Takaharal[T9M iKataoka et alJ[T999h . 
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Although the simplest model successfully explains the vari- 
ability and the broadband spectra from radio to TeV 7 ray 
band to some extent, recent improvements of atmospheric 
Cherenkov TeV 7-ray telescopes such as H.E.S.S, MAGIC 
and VERITAS make it possible to provide high quality spec- 
tra in relatively short exposures (see, e.g., Baixeras 200^ 
iHinton 2004; Weekes et al.ll2002l) . Indeed, the recent obser- 
vations with those telescopes imply the conventional model 
may no longer be applicable. In particular, the "orphan" 
flares detected from IBS 1959-H650 (Kra wczynski et alJ2004 
and Mrk 421 (Blazejowski et al. 2005), during which the 
flux of TeV 7-ray flared up without counterparts in X-ray 
band, are difficult to explain within the context of simple 
one-zone SSC models. Also, the detection of short variabil- 
ity time scale of TeV 7-ray emission from Mrk 421 ('^ 10 
minutes; Gaidos et al. 1996) and PKS 2155-304 a few 
minutes; Ahar onian et al.l 120071) may require unprecedent- 
edly l arge Doppler factor (5 > 50 (e.g.. lGhisellini & Tavecchio 
I2OO8I) . which are much larger than ones consistently estimated 
from multi-wave length analysis (e.g.. iTavecchio et al.lll998l; 
iKubo et al.l[T998h . 

In order to understand the particle acceleration and the non- 
thermal emission from TeV blazars, more sensitive observa- 
tions are needed. In particular, hard X-ray observations above 
10 keV are important because the emission in this energy 
band reflects the behavior of the most energetic electrons and 
is a sensitive probe of their acceleration, cooling and escape 
rates. In this paper, we present the results of the observa- 
tions of Mrk 421 with Suzaku, which is one of the most suit- 
able observatories for exploring such complex behaviors with 
its high sensitivity and wide-band coverage from the soft X- 
ray to hard X-ray band. The observation was performed in 
2006 April 28-29, at a time when the object was undergoing 
an outburst and showed a high level of activity. In § 2 and 
§ 3, we describe our Suzaku observations and the data reduc- 
tion procedures. Analysis and results are shown in § 4, and 
are discussed in the following § 5, with summary in § 6. In 
the present work, we use the data products from the Suzaku 
pipeline processing version 2.0. The data reduction and anal- 
ysis are done using HEADAS 6.3.1 and the spectral fitting 
is performed with XSPEC 11.3.2. Regarding the notation 
for both photon and electron energy, e and ly are for photon 
energy and frequency, and E and 7 are particle energy and its 
Lorentz factor, respectively. 

2. OBSERVATIONS 



The X-ray Observatory Suzaku jMitsuda et al.ll2007l) . de- 
veloped jointly by Japan and the US, has a scientific payload 
consisting of two kinds of co-aligne d instruments; the XIS 
jKovamaetal. 2007) and the HXD dTakahashi et al.1 120071: 
iKokubun et al.ii2007i) . The XIS consists of four X-ray sen- 
sitive CCD cam eras which ar e located in the foci of X-ray 
telescopes (XRT; Serlemitsos et al.l2007l) . The HXD is a non- 
imaging detector system which covers the hard X-ray band- 
pass of 10-600 keV with PIN sihcon diodes (10-60 keV) and 
GSO scintillators (40-600 keV). 

We performed the Suzaku observation of Mrk 421 starting 
from 2006 April 28 (MJD 53853) 06:46 UT through April 
29 (MJD 53854) 06:30 UT Figure [U shows a long term fight 
curve of Mrk 42 1 obtained with RXTE ASM (2-10 keV, 
iBradt et aTl Il993h . which includes the time interval of our 
Suzaku observation. As shown in the figure, Mrk 421 ex- 
hibited high X-ray flux with a peak value reaching ~ 2.5 
counts/sec around MJD = 53853, which is corresponding to 
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Fig. 1.— A long-term light curve of RXTE ASM (2-10 keV) around the 
Suzaku observation. 

~ 30 mCrab. All the XIS sensors were operated with 1/4 
window option in order to reduce possible pile-up effects. 
In this mode, the restricted detector regions corresponding to 
17'. 8 X 4'. 5 in the sky are read out every 2 sec. 

3. DATA REDUCTION 

For the XIS analysis, we retrieve "cleaned event files", 
which are screened with standard event selections. We fur- 
ther screen the events with the following criteria: (1) cutoff 
rigidity larger than 6 GV/c and (2) elevation angle larger than 
5° and 20° from the Earth rim during night and day, respec- 
tively. We do not use XIS 1 data in the following analysis since 
it suffered from telemetry saturation for almost the entire ob- 
servation period. 

The XIS events are extracted from a circular region with a 
radius of 176 pixels 3') centered on the image peak. This 
extraction circle is larger than the window size and therefore 
the effective extraction region is the intersection of the win- 
dow and this circle. Since the background is entirely negli- 
gible throughout the observations (less than 1 % of the sig- 
nal even at 9 keV), we do not subtract background from the 
data. The source was so bright during these observations that 
the XIS suffered from photon pileup at the image center even 
with the 1/4 window option. In the present study, we exclude 
a circular region of the radius of 20 pixels at the image cen- 
ter from the event extraction region to minimize these effects. 
After excluding the central region, the systematic error due to 
the pile-up effects included in the flux is estimated to be less 
than 2 %. 

In order to take into account the complex shape of the 
event extraction region (an annulus truncated with the window 
size), we calculate the response matrices (RMF) and the ef- 
fective area (ARE) for each XIS sensor by using xi s rmf gen 
and xissimarfgen jlshisaki et al. 2007), respectively. The 
xissimarf gen is based on Monte Carlo simulation and is 
designed to properly handle the geometrical shape of the ex- 
tracted region. Since the three EI CCDs have almost identical 
performance, we sum their data to increase statistics in the 
spectra. Accordingly, we calculate the corresponding RME 
and ARE response/area files by adding those for XIS. 

Eor the HXD data, "uncleaned event files" are screened 
with standard event screening criteria. We exclude events dur- 
ing (1) SAA passages and (2) Earth occultation, and also (3) 
those with cutoff rigidity less than 6 GV/c. To estimate accu- 
rately the Non-X-ray Background (NXB) is the key to realize 
the high sensitivity observation for HXD. In this analysis, we 
use bgd_d (for rev 1 .2) for the NXB of HXD/PIN because the 
current public background model (bgd_a) is found to overesti- 
mate the PIN background in the period during 2006 Mar 27 to 
May 13 tEukazawa et al...2009.) . Note that we have to correct 
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Fig. 2. — Light curve of HXD/PIN during Earth occultation. In the upper 
panel, the obtained data and the background model in the energy band of 
15-60 keV are compared. The bottom panel shows the differences between 
them. 



the effect of dead time by hxddtcor for this bgd_d as well 
as the source event files. T he current NXB model is shown 
to be accurate within ^3 % dFukazawa et al.l l2009). Since the 
HXD/GSO detected no significant signals above 2 a assum- 
ing 5 % systematics for the NXB estimation, below we report 
results obtained from the XIS and the HXD/PIN. 

The accuracy of the background model is evaluated by com- 
paring the count rate during the Earth occultation between the 
observation data and the NXB model, since the NXB model 
is constructed based on the Earth occultation, which is de- 
fined by the elevation angle from the Earth rim less than -5° 
(ELV < -5). In Figure |2] we compare the light curves of 
HXD/PIN data in 15-60 keV band actually taken during the 
Earth occultation and of the NXB model. Each data point 
is about 1 ksec accumulation. Figure |2] shows that the NXB 
model is well reproduced with an accuracy of less than 5 % 
of the NXB even for 1 ksec observation. Therefore, the un- 
certainty due to the misestimation of the NXB does not affect 
our results, because the flux is more than 30 % of the NXB 
above 25 keV even in the lowest flux period throughout the 
observation, as we will see below. 

Since the NXB model does not include contributions from 
the CXB (Cosmic X-ray Background), a simulated spec- 
trum of the CXB is added to the NXB model. Following 
iGruber et al.l (119991) . who reanalyzed HEAO-1 observations 
in the 1970s, our spectral model for the extragalactic back- 
ground radiation is chosen as 



~r ~ keV exp 

at 



CkeV , , - 

ph s 

En 



keV 



-I -2 - 

cm sr 



(1) 



where e^ev = e/1 keV and ep = 41.13. We estimate the 
expected CXB signal in the HXD/PIN spectrum using 
the latest response matrix for spatially uniform emission, 
ae_hxd_pinf latel_20070914 . rsp. The contribu- 
tion from the CXB flux is estimated to be 5 % of the NXB, 
which is comparable to the current systematic uncertainties of 
the NXB model itself. 
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Fig. 3.— The sum of data from three XIS sensors and HXD/PIN light 
curves in the energy ranges of 0.5-2.0 keV, 2.0-10.0 keV, 12-20 keV, 20- 
30 keV and 30-60 keV. The dot-dash line coiTesponds the 10% of NXB for 
30-60 keV band. 

4. ANALYSIS 
4.1. Li^M Curves 

The time history of Suzaku observation obtained with XIS 
and HXD is shown in Figure[3]for five energy bands. The flux 
in each energy band decreases by factor 2-A during the first 
half of the observation and then turns into an increase starting 
at around the middle of the observation. Thanks to low back- 
ground characteristics of HXD, the variability with time scale 
of as short as 20 ksec is clearly detected even at the highest 
band, 30-60 keV. The dot-dash line drawn in the bottom of 
the figure represents the 10 % of the estimated NXB for the 
highest energy band, whereas the contribution from CXB in 
the same energy band is about 0.002 counts/sec. This also 
supports the significant detection of the signal from Mrk 421 
up to 60 keV in such short time intervals. 

The change in flux becomes larger in higher energy bands, 
which is consistent w ith previous studies of Mrk 421 (e.g., 
iTakahashi et al.llT996l 1200 0). The count rate decreases from 
100 count s"^ down to 70 count s"' in the lowest energy 
band (0.5-2.0 keV in the XIS data), while it drops from 0.07 
counts/sec down to 0.02 counts/sec in the highest energy band 
(30-60 keV in the HXD data), that is, the higher energy band 
is more than twice as variable as the lower energy band, even 
considering the systematic and statistical errors. 
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4.2. Suzaku Spectral Analysis 

In order to study the time evolution exhibited in the light 
curves more quantitatively, we divide the XIS and HXD 
data into 14 time intervals and perform spectral fitting to 
XIS+HXD spectra derived from each time interval. Raw 
HXD/PIN spectra during the highest and the lowest flux pe- 
riods are shown in Figure HI During the observation, the flux 
below 20 keV is always higher than that of the NXB. At the 
lowest flux period (lower panel of Figure|4|, flux level is com- 
parable to the NXB ai-ound 25 keV and about 30 % of the 
NXB level above 25 keV. Since the syst ematic error of the 
NXB is about 5 % jFukazawa et al.lf2009l) . which is six times 
less than the flux level measured at the lowest flux period, 
the results of present analysis are not affected by the sys- 
tematic uncertainties of the NXB. In the following spectral 
analysis, the column den sity is fixed to the Galacti c value: 
Nh = 0.0161 X 10^2 cm-2 (iLockman & SavagellT995h and the 
cross normalization between XIS and HXD/PIN is also fixed 
to 1.13 (Kokubun et al. 2007). 

Wide band spectrum obtained by Suzaku enables us to mea- 
sure time-resolved behavior of the continuum spectrum, in- 
cluding even small deviations from a simple power law model. 
Firstly, in order to see how the spectra deviate from power law 
function, we restrict the fitting range to 2-5 keV. The results of 
fitting are shown in Figure|5] and as is clearly seen in the ratio, 
the spectra gradually curve downwards for both cases, mean- 
ing the power law index increases with increasing energy. 
Such a departure from a simple power law model appears to 
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Fig. 4. — The raw spectra of HXD/PIN at the brightest period (upper panel) 
and at the lowest period (lower panel). Only 750 sec accumulation time is 
used for the brightest period and 1.7 ksec for the lowest period. The open 
circles, open squares and filled circles show the raw data, background (NXB 
+ CXB), and background-subtracted data, respectively. 
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Fig. 5. — XIS+PIN spectra in the energy range 0.4-60 keV are shown with 
the best-fit models obtained by fitting the XIS spectrum (2-5 keV) alone with 
power law model. The data points below/above 10 keV correspond to XIS 
and PIN data respectively. The lower panels show the ratio between the data 
points and the model values for the highest and the lowest periods. 



be a general feature of X-ray spectra of BL-Lac type blazars 
with the peak of the s ynchrotron compone nt in the UV to soft 
X-ray band (see, e.g.. lPerlman et al]l2005l) . It is apparent that 
the spectrum changes with time, even in this restricted energy 
range: the spectral indices of power law model (2-5 keV) 
significantly change from 2.05±0.02 to 2. 27 ±0.03, with the 
flux-index correlation in the sense where the soft X-ray spec- 
tra soften as the flux decrea se. This is consistent with the 
results published so far ('e.g.. lTakahashi et"anil996h . Clearly, 
sensitive Suzaku observations indicate spectral variability on 
short time scales (a few hour), and provide an opportunity to 
derive important constraints in the context of more physical 
"synchrotron model," where the data are directly fitted to a 
particle spectrum radiating via the synchrotron process. This 
treatment is covered in the following sections. 

4.3. Interpreting the X-ray Spectrum in the Context of 
Physical Synchrotron Model 

We can obtain constraints on the physical parameters of 
the synchrotron emission, such as B and 7max, via fitting the 
spectra with the recent numerical model (Tanaka et al. 2008). 
Hereafter we simply call this model "synchrotron model." In 
the model, we assume an electron spectrum has a form of 
power law with exponential-type cutoff at the maximum en- 
ergy Emax (= Tmax'WeC^), namely. 



dNe 

—— = NeoE ■ exp 
dE 



(2) 



Here, s is the energy index of electron distribution, while f3 
and Neo (electrons cm""* eV~') determine the shape of high 
energy end of the pop ulation and the normalization, respec- 
tively. As described in Kataok a et al.l ([199^, we assume the 
spherical and uniform emitting region. The total power per 
unit frequency P(e) can be written as 



Aot(e) ^ I F{- 



dE 



dE. 



(3) 
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Here e and E parameters indicate the photon and electron en- 
ergy, respectively, and the function F{x) is defined as 



F(x)=x I KspiOd^ 



(4) 



where K^/iiO is the modified Bessel function of 5/3 order. 
When pitch angles are isotropic, the characteristic photon en- 
ergy, is given as 

2 



:5.43 



B 



0.10 G 



1 TeV 



keV. 



(5) 



This model has three free parameters; s, H = E^.^^B^I^ [GeV 
G'/^] and the normalization of a spectrum. We note here that 
this synchrotron model is in fact based on the distribution 
of energies of radiating particles that is essentially a power 
law with an exponential (or hyper-exponential) cutoff. Such 
a distribution is similar to that which has been calculated to 
stochastic acceleration via the Fermi 11 process - such as, 
for instance, via scattering in m agnetic turbulent regions in 
plasma (see, e.g. JStawarz & Petrosian 2008). 

In the case of emission from a jet, we need to consider 
the effect of relativistic beaming, since the fitted value of 
n (Hobs) is not given in the jet's frame (Iljet), but instead 
Hobs = Iljet X We adopt the Doppler (beaming) factor 

(5=10: this value was directly measur ed by VLBl obser- 
vations (e.g.. lVermeulen & Coheniri994l) which is consistent 
with the result obtained by simultaneous multi-band analysis 
from radio to 7-ray (QC ubo et al. 1998). We note that more re- 
cent work by L ichti et al.i (2008) which analyzes the big out- 
burst followed by some weeks with respect to our Suzaku ob- 
servation and suggests 5= 15 by modeling the SED. However, 
these differences are not the scope of our paper and will not 
affect our discussion. 

When we perform the spectral fitting for each data set with 
this model with (3 fixed to unity, almost all the individual spec- 
tra are well fitted by this synchrotron model. In the fitting, 
we exclude the bandpass 1.5-2.5 keV and 8.0-10.0 keV be- 
cause there remains some calibration uncertainties of the XIS 
detectors. Moreover, we fix the energy index of electron dis- 
tribution (5) to be 2.0: then, the reduced chi-square (x^) is 
ranging from 0.86 to 1.57 with 57 degrees of freedom (dof). 




2.0 5.0 10.0 20.0 

Energy [keV] 



Fig. 6. — The same as Figure|5]but fitting with simple synchrotron model 
in entire energy band 0.4—60 keV. 
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Fig. 7. — The results of fitting time-resolved spectra with simple syn- 
chrotron model. The top panel represents the evolution of Elobs values in 
units of GeV G'^^ obtained by fitting each spectrum in the case of the energy 
index of electron spectrum fixed to s = 2.0 (red), s = 2.5 (green) and .v = free 
(black). The error bars include only statistical Icr errors and those for green 
and red diamonds are too small to be seen. The middle panel shows the his- 
tory of the energy index of electron spectrum in the case of i = free. In the 
bottom panel, we show the xi values corresponding to 59 and 60 d.o.f. 

In some time intervals, the xi values are substantially larger 
than in others and the model will be rejected. However, rel- 
atively large values of xi are partially caused by the insuffi- 
cient calibration of the XIS/XRT rather than the inappropriate 
modeling of the spectra. Figure |6] shows the typical examples 
of the results for the highest- and the lowest-flux periods. For 
those, we obtain Hobs = 243 ± 2 (xl = 1 .42 for 57 dof) and 
149 ± 1 (x?, = 1 . 1 8 for 57 dof) [GeV G'/^] for the highest and 
the lowest period, respectively. 

We summarize the fitting results of each spectrum by the 
"synchrotron" model in Figure [T] The top panel of Figure [T] 
shows the evolution of Hobs as a function of time, where the 
data points in red color correspond to the case of s = 2.0. 
With s = 2.0, the Hobs = £'maxB'''^(5'/^ varies largely as the flux 
changes. In order to study how the energy index of electrons 
affects the results, we also fit those spectra with s = 2.5, which 
is shown by the green markers in the Figure [T] The Hobs val- 
ues in the case of s = 2.5 are significantly different (higher) 
from those in the case of i = 2.0. In both cases. Hobs param- 
eters change by a factor of 2. On the other hand, when we 
treat i as a free parameter. Hobs does not change significantly, 
while s increases from 1.95 to 2.55 in the decreasing phase. 
Statistically, xi for the case of s as free is the lowest in any 
time intervals compared with in the case of fixed s by only 
an additional degree of freedom, except for the first and the 
second time intervals. The interpretation of those results is 
discussed in § 15.11 

The spectra are fitted fairly well with the model based on 
synchrotron emission. In Figure [8] we plot the spectra in 
v-vFi, space for all the data sets. All the spectra have con- 
vex shape. The shape of these spectra change significantly. 
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Fig. 8. — SED for each good time interval fitted by "synchrotron" model. 



and that is quite pronounced in the first interval (#1). In the 
first half of the observation, the flux decreases as the syn- 
chrotron peak frequency gradually shifts to lower band, that 
is, to ^ 5 keV at the highest period (#1) and below 0.4 keV 
at the lowest period (# 9). On the other hand, in the second 
half, the synchrotron peak frequency does not seem to change 
apparently, as the flux increases. 

Since X-ray spectra of blazars are often fitted to a spectral 
model involving a broken power law, for comparison, we at- 
tempted such a model. In all cases, the broken power law 
(involving three free parameters plus normalization, namely 
index before the break, after the break, and the break energy) 
returns a worse fit than the synchrotron model (where we used 
i fixed at 2): for instance, for period # 1, x?. is 1.42 for 57 dof 
for the "synchrotron" model while it is 2.28 for the broken 
power law; for period #2, the values are respectively 1.27 vs. 
1 .40: the broken power law model appears to have a break that 
is too sharp to fit the data. As it will become important in the 
discussion below, most relevant is the comparison of the two 
models for the period # 9: here, the values are respectively 
1.18 vs. 1.90. 

4.4. Spectral Analysis of the Steady and the Variable 
Components 

One possible approach towards the exploration of the spec- 
tral evolution during the observation is to separate the X-ray 
data into a "steady" and a "variable" components. In reality, 
"steady" means slowly-variable, variable on time scales much 
longer than typical X-ray observations lasting for one or sev- 
eral days, since historically the source flux does drop well 
below the level of our "steady" level. We associate the "vari- 
able" component with a newly injected electron distribution 
which drives the flare activity. We consider the "steady" com- 



ponent to correspond to the lowest state spectrum (# 9), which 
is illustrated in Figure|4]and well described by the synchrotron 
model in Figure|6] and the "variable" component to be the dif- 
ference of individual spectra obtained at other epochs minus 
that of the lowest state spectrum. 

In order to unveil the characteristics of "variable" compo- 
nent, we fit each of the subtracted spectra with simple power 
law model. The photon index of the subtracted spectrum is 
found to be harder, F ^ 1.8, but it also shows a large discrep- 
ancy between data and model and therefore is not well fitted 



5.0 10.0 
Energy [keV] 



Fig. 9. — Power law and broken power law fit for a subtracted spectrum. 
This is a case of second time interval (# 2) out of 14. The middle and bottom 
panels represent the ratio between data and model for simple power law and 
for broken power law model. 
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TABLE 1 

Fit results of the "variable" components with a broken 
power law function " 



Cbrk 



2^10 keV flux'' xlidof) 



1 

2 
3 
4 
5 
6 
7 

11 
12 
13 
14 



1 63+""2 

1 70+0-02 
1 63+003 
1.68: 
1.80 



i+S:l 



1.17: 
1.46: 
1.45: 

1.58 



1^ 

'-0.05 



3 6+<' 3 

•'"'-0.3 
n 4+0.3 

-0.2 

3.4+0-2 

^•^-0.5 

2 9+fl.s 



J 9+0.5 

2.0+^-^ 
3 0+^-^ 

^■^-0.3 



2 jj+0.04 

2 1 6+^-^^ 
9 , 0+0.04 
'^■^-'-0.04 

7 nq+004 
-!!■(•!( 

T 91 +0.10 
^■^ -0.08 

9 o9+<).20 
■^■-^ -0.16 

1 Q7+fl.I4 

9 m+8:M 

'^■'-"-0.07 



4.12 X 10- 


10 


1.25 (54) 


2.49 X 10" 


10 


1.08 (54) 


2.16 X 10" 


10 


0.87 (54) 


2.10 X 10" 


10 


1.23 (54) 


1.50 X ur 


10 


0.97 (27) 


7.90 X 10" 


11 


0.81 (27) 


4.88 X 10" 


11 


1.17(27) 


3.75 X nr 


11 


0.72 (27) 


4.70 X ur 


11 


1.14(27) 


6.91 X ur 


11 


0.75 (54) 


1.21 X 10" 


10 


0.80 (54) 



^ The evolution of these parameteis i,s shown in Figure 1 1 01 

^ The Flo and Phi indicate the photon indices below and above the break energy ebii< 
The break energy €brk is in units of keV. 
The 2-10 keV flux is in units of erg s"' cm"". 



by this model. Next we attempt a broken power law model 
and find that this model is acceptable for all the time intervals. 
An example of the fitting is shown in Figure|9]comparing with 
that of simple power law model, corresponding to the epoch 
# 2, and xi is significantly improved from 3.13 (56 dof) to 
1.08 (54 dof). Subsequently, we fit all of the intervals with 
such a broken power law spectrum : Figure[TO]shows the time 
history of parameters, that is, the photon index of low energy 
band (Fio) and of high energy band (Fhi), and the break energy 
(cbrk), which are summarized in Table [T] The vicinity of the 
epoch # 9 are excluded because the relatively large statistical 
errors of energy bins make it difficult to extract meaningful 
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Fig. 10. — Evolution of "variable" spectra. Top and middle panels show 
the history of the photon indices for lower (black; Fio) and higher (red; Fhi) 
energy band and the break energy (ebin) in units of keV. We consider the 
"steady" component to coirespond to the lowest-flux state (# 9). The vicinity 
of # 9, that is, # 8 and # 10 are excluded because of relatively poor statistics. 



Bottom panel shows that of the 2-10 keV flux in units of erg s cm 



-2 



values. The 2-10 keV flux of the "variable" component de- 
creases from 4.1 X 10"'° to 4.9 x 10"" erg s"'cm"^ (and of 
course to zero for interval # 9, by construction) with the e- 
folding decaying time of 1.7 x 10"* sec, that is, less than 5 
hours. 

These plots lead to two important results : (1) the spectral 
shape of "variable" component does not vary across entire 
observation and only the normalization changes and (2) the 
photon indices Fio and Fhi are about 1.6 and 2.1 respectively, 
that is, the differences between them are equal to about 0.5, 
with the break energy in photon space remaining relatively 
constant, at ebrk — 3 k eV. All this has important implications 
which we discuss in § 15.21 

It is interesting to compare the lowest spectrum - that corre- 
sponding to the # 9 epoch - to the remaining spectra. This 
is essentially a comparison of the "variable" and "steady" 
components. In Figure [TT] we plot both spectra at the same 
graph in the v-vFi, space; there, we illustrate the highest flux 
period (# 1 minus # 9) and the end of decaying period (# 
6 minus # 9) for the "variable" components. The spectral 
shapes of the "steady" and "variable" components are quite 
different, specifically the spectral indices of "variable" com- 
ponents, both lower and higher energy band, are harder than 
these of "steady" component. Perhaps the most striking is the 
fact that the spectral shapes of "variable" components do not 
vary significantly, and only the normalization changes. 

5. DISCUSSION 

In the previous section, we performed detailed spectral 
analysis of TeV blazar Mrk 421 obtained by Suzaku in the 
period of the large flare in 2006 April 28-29 and found that 
there are two different interpretations for their variability. Be- 
low, we consider the physical scenarios that might correspond 
to the two phenomenological models above. 

5.1. Spectral Variability due to the Change o/IIobs 

In § 14.31 we successfully fitted the evolution of the X-ray 
spectrum by the synchrotron model. In the case where the en- 
ergy index of electron population is allowed to be free, Ilobs 
parameter does not change significantly throughout the obser- 
vation and the change of electron energy index s appears to 
be the factor responsible for the spectral change. The electron 
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Fig. 1 1. — Simultaneous plots of "variable" and "steady" components for the highest flux period (epoch # 1, left panel) and the end of decay period (epoch 
# 6, right panel). The red and black data set are for the "variable" components and the total (="steady"+"variable") spectra. The blue data set con'esponds to the 
lowest-flux peiiod and is common to two plots. 



energy index itself is unlikely to vary significantly, so if the 
spectral evolution is due to the change of single electron pop- 
ulation, the observed variability must be due to the change of 

Hobs- 

Using equation (2.31) in llnoue & Takaharal (119961) . Hobs 
can be written as 



Hobs = B'''^(5'''^7tnax'«eC^ = V^BnieC 



9eS 



80(MB+Msync)a'TC 



1/2 



(6) 

where v.,, ub, Usync, ctj and ^ are the shock velocity in the ob- 
server's frame, the energy density of magnetic field and syn- 
chrotron seed photons, Thomson cross section, and the so- 
called gyro-factor (the ratio of the mean free path for scatter- 
ing with magnetic disturbances compared to the Larmor ra- 
dius), respectively. This equation assumes that the maximum 
energy of electron distribution is determined by equating the 
radiative cooling time to the accelerating time. From this 
equation, the variability of Hobs can be explained by changes 
of one or more parameters such as Vj, B (thus also ub), Msync, 
^, or 5. The degeneracy between the parameters can be bro- 
ken by comparing the temporal and spectral behavior of the 
X-ray time series against the TeV flux history. For instance, 
if the magnetic field B is mainly responsible for the changes 
of Hobs, it would require for the ratio of the synchrotron lu- 
minosity to the inverse Compton luminosity, L^yac/Lic, sig- 
nificantly change as a function of time. However, such multi- 
wavelength studies require simultaneous X-ray and TeV data, 
which are not currently available for this Suzaku observation. 

5.2. Spectral Variability due to Superposition of Steady and 
Variable Components 

As an alternative and preferred scenario to explain the vari- 
able behavior of Mrk 42 1 , we introduce two separate electron 
populations, one is "steady" and the other is "variable" corre- 
sponding to the steady and variable componsnts of the X-ray 
spectra. There, a fresh electron distribution, distinct than the 
pre-existing one, is newly injected into the same (or a dif- 



ferent) region of the jet and is responsible for the observed 

variable portion of source flux. 

One of the important results in § 14.41 is that the photon 
indices Fio and Fy are measured to be '-^1.6 and ~ 2.1 
for "variable" components. Surprisingly, the result suggests 
that electrons form the power law distribution of energy in- 
dex ^ 2, which is predicted from the standard shock accel- 
eration theory for both non - relativistic and relativistic case 
(iBlandford & Ostrikeii[T97l iKirk et all 12000) and then suf- 
fer from synchrotron cooling. As for the variability, we have 
shown that the spectral shape of "variable" components does 
not vary throughout the observation and the break energy ap- 
pears to stay constant around ebik ^ 3 keV. While it is pos- 
sible to adjust other parameters simultaneously to keep ebrk 
constant, this would require simultaneous changes of several 
parameters, amounting to "fine-tuning" of parameters. The 
simplest interpretation of this observational result is that the 
variability of the X-ray flux is only due to the change of the 
injection rate of electrons, A^^) in Eq. (|2]i. Since the break en- 
ergy in the synchrotron spectrum (ebrk) is constant, B, 5 and 
7bik are also expected to stay constant. With this, the syn- 
chrotron luminosity is proportional only to the energy density 
of the particles: 



Lsync OC U^UbS'* CX Me(5^7bit OC A^^ 



eO- 



(7) 



We thus decompose the X-ray spectrum of Mrk 421 into 
a "steady" component described as due to an exponentially 
cutoff power law distribution of electrons, and the "variable" 
component, well-described as a broken power law photon 
spectrum. Here we can make a suggestion of physical pic- 
ture for these two components: the "variable" component is 
due to electrons efficiently accelerated in relatively small re- 
gions, such as localized shocks, via Fermi I process. In ad- 
dition to that, a relatively "steady" component contributes as 
well: one possibility would be that the observed spectrum of 
this component is produced by shocks located at larger dis- 
tance along the jet from the black hole. Such picture, invoking 
an internal shock scenario, was suggested by .Tanihata et akl 
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( 12003b . There, the resolved, rapid flares are due to collisions 
of pairs of shells at a characteristic distance Dch from central 
engine, while those colliding at larger distance than D^i make 
up the underlying, slowly variable component. Alternatively, 
the "steady" component might be associated with another pro- 
cess: the electrons escaping from the shock region into more 
extended volume are re-accelerated via Fermi 11 process on 
turbulence sites, as would be in the scen ario described by 
several auth ors ( Virtanen & Vainio 2005^ iKatarzvnski et alj 
l2006t IStawarz & Petrosianii200 8). This process operates in 
a larger volume, and thus the longer variability time scale. 
This scenario is particularly attractive as it is strongly sup- 
ported by the spectral shapes and relative variability time 
scales of the two putative components. Specifically, the pho- 
ton sgectrum calculated for the diffuse acceleration process 
by IStawarz & Petrosian ( 2008) is similar to that derived via 
the synchrotron model, described and applied by us above. 

With regard to the "variable" component, we infer an im- 
portant implication of the fact that the observed break energy 
seems not to change throughout the observation. There, the 
break seen in the photon spectrum might be due to a competi- 
tion between acceleration and cooling of relativistic electrons 
and expressed as 



3otpC 



7brk ■■ 



4(MB+Msync)0'T^ 



(8) 



where R is the size of emitting region (llnoue & Takaharal 
11996^. To be exact, Msync in Eq.dHJ is actually an increasing 
function of A^^-o via Eq. (Ql and this implies Cbrk should increase 
as the flux of "variable" component decreases. However, the 
fact that ebrk does not seem to change throughout the decaying 
phase of "variable" component suggests that the energy den- 
sity of magnetic field is larger than that of soft seed photons 
{ub ^ Msync) for "variable" components. This in turn suggests 
the synchrotron loss is dominant over the inverse Compton 



^ 1 ksec. 

Our detailed, time-resolved spectral analysis leads to two 
different scenarios describing the variability of Mrk 42 1 . One 
is the conventional, one-component picture: in order to in- 
terpret the spectrum, we introduce the synchrotron model for 
the fitting function and find that the variability is due to the 
change of Hobs, the products of Z?'/^, (5'/^ and 7maxWieC^. The 
variability of Hobs can be explained by changes of one or more 
of them. However, since we cannot break the degeneracy only 
this observation, simultaneous X-ray and TeV 7-ray observa- 
tion are required. The other scenario invokes also a second, 
separate electron distribution that is responsible for the vari- 
ability and which has energy index 5^2; this distribution is 
distinct than the one responsible for the steady component. 
Here, the rapidly variable component might be due to local- 
ized shock (Fermi 1) acceleration, while the more steady com- 
ponent might be due to the superposition of shocks located at 
larger distance along the jet or due to a larger-scale, stochas- 
tic (Fermi 11) acceleration on turbulence sites in the shocked 
plasma. 

It is not possible to distinguish which scenario correctly de- 
scribes the emission mechanism solely via this observation. 
However, the multi-component scenario provides a meaning- 
ful hint to disentangle the puzzle of the emission mechanism, 
such as the "orphan" flares and the very short time scale 
minutes) of variability observed in TeV 7-ray band. In or- 
der to investigate the emission mechanism of TeV blazars fur- 
ther, simultaneous broad-band, high-sensitivity observations 
are needed, especially in X-ray and GeV/TeV 7-ray band fa- 
cilitating a comparison of the "variable" components over the 
entire broad-band blazar spectra. Those show a promise to- 
wards a progress for understanding the acceleration mecha- 
nism in the relativistic jet. Such observations will be possi- 
ble using Suzaku, Fermi GeV 7-ray telescope and advanced 
Cherenkov 7-ray telescopes. 



loss (Ub/Us 



sync " 



-^sync,var 
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6. SUMMARY 

We performed a detailed analysis of bright TeV blazar 
Mrk 421 data collected with Suzaku during the large flare in 
2006 April 28-29. During the observation, the flux varied 
from 3.7 x 10"'° to 7.8 x 10""' erg s"'cm"2. Thanks to the 
high sensitivity of Suzaku, we obtain the time-resolved spec- 
trum of wide energy band 0.4-60 keV in intervals as short as 
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